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ABSTRACT 

The transition redshift (deceleration/acceleration) is discussed by expanding 
the deceleration parameter to first order around its present value. A detailed 
study is carried out by considering two different parameterizations: q = qo + 
qiz and q = qo + qiz{l + z)~^, and the associated free parameters {qo,qi) are 
constrained by 3 different supernova samples. The previous analysis by Riess 
et al. [ApJ 607, 665, 2004] using the first expansion is slightly improved and 
confirmed in light of their recent data {Gold07 sample). However, by fitting the 
model with the Supernova Legacy Survey (SNLS) type la sample we find that the 
best fit to the redshift transition is Zt = 0.61 instead of Zt = 0.46 as derived by the 
High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample 
is also in good agreement with the Davis et al. sample, Zf = 0.60^q;^^ (1<^)- Such 
results are in line with some independent analyzes and accommodates more easily 
the concordance fiat model (ACDM). For both parameterizations, the three SNe 
type la samples considered favor recent acceleration and past deceleration with 
a high degree of statistical confidence level. All the kinematic results presented 
here depend neither on the validity of general relativity nor the matter-energy 
contents of the Universe. 

Subject headings: Cosmology: observations - decelerating parameter - kinematic 
model 



1. INTRODUCTION 



It is now widely believed that the universe at redshifts smaller than unity underwent 
a "dynamic phase transition" from decelerating to accelerating expansion which has been 
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corroborated by several independent analyzes. In the context of the general relativity theory 
such a phenomenon can be interpreted as a dynamic influence of some sort of dark energy 
whose main effect is to change the sign of the universal decelerating parameter q{z). 

The most direct observation supporting the present accelerating stage of the Universe 
comes from the luminosity distance versus redshift relation measurements using supernovae 

(SNe) type la (Riess et al. 1998; Perlmutter et al. 1999; Riess et al. 2004, 2007; Astier et 
al. 2006; Wood-Vasey et al 2007; Davis et al. 2007; Kowalski M., et al. 2008) initially 
interpreted in light of ACDM scenarios using either background or inhomogeneous luminosity 
distances (Santos et al. 2008; Santos & Lima 2008). However, independent theoretical 
observational analyzes points to more general models whose basic ingredient is a negative- 
pressure dark energy component (for review see, Padmanabhan 2003; Peebles & Ratra 2003; 
Lima 2004; Copeland et al. 2006). 

The convergence of many high-quality experimental results are now strongly suggesting 
that the cosmic picture at present is the following: the spatial geometry (curvature) is flat, or 
more approximately flat, and the dynamics is governed by a component called dark energy, 
3/4 of composition, and 1/4 for matter component (baryons plus dark). Among a number 
of possibilities to describe this dark energy component, the simplest and most theoretically 
appealing way is by means of a positive cosmological constant A. Others possible candidates 
are: a vacuum decaying energy density, or a time varying A-term (Ozer & Taha 1987; Freese 
et al. 1987; Carvalho et al. 1992; Lima & Maia 1993; 1994; Maia & Lima 1999; 2002; Lima 
1996; Ovcrduin & Coopcrstock 1998; Cunha et al. 2002a,b; Cunha & Santos 2004; Alcaniz & 
Lima 2005; Carneiro & Lima 2005; Fabris et al. 2007), a time varying relic scalar field slowly 
roUing down its potential (Ratra & Peebles 1988; Prieman et al. 1995; Caldwell et al. 1998; 
Saini et al. 2000; CaldweU 2000; Carvalho et al. 2006), the so-called "X-matter", an extra 
component simply characterized by an equation of state Px = <^Px (Turner & White 1997; 
Chiba et al. 1997; Alcaniz & Lima 1999, 2001; Kujat et al. 2002; Alcaniz et al. 2003), the 
Chaplygin gas whose equation of state is given hy p = —A/p where A is a positive constant 
(Kamenshchik et al. 2001; Bcnto et al. 2002; Dev et al. 2003; Bento et al. 2003; Cunha 
et al. 2004; Lima et al. 2006a,b), among others (Lima & Alcaniz 1999; Chimento et al. 
2001; Freese & Lewis 2002; Pavon & Zimdahl 2005). For SFC and XCDM scenarios, the uu 
parameter may be a function of the redshift (see, for example, Efstathiou 1999; Cunha et al. 
2007), or still, as it has been recently discussed, it may violate the dominant energy condition 
and assume values < — 1 when the extra component is named phantom cosmology (Caldwell 
2002; Lima et al. 2003; Perivolaropoulos 2005; Gonzalez-Diaz & Siguenza 2004; Lima & 
Alcaniz 2004; Santos & Lima 2008). It should be stressed, however, that all these models 
are based on the validity of general relativity or some of its scalar-tensorial generalizations. 
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On the other hand, Turner & Riess (2002) have discussed an alternative route - some- 
times called kinematic approach - in order to obtain information about the beginning of 
the present accelerating stage of the Universe with no assumption concerning the validity of 
general relativity or even of any particular metric gravitational theory (in this connection sec 
also Weinberg 1972). Although considering that such a method does not shed light on the 
physical or geometrical properties of the new energetic component causing the acceleration 
it allows one to assess the direct empirical evidence for the transition from decelerating to 
accelerating in the past as provided by SNe type la measurements. In their preliminary anal- 
ysis it was found that the SNe data favor recent acceleration {z < 0.5) and past deceleration 



More recently, with basis on the same approach, the High-z Supernova Search (HZSNS) 
team have obtained Zt — 0.46 ± 0.13 at la c.l. (Riess et al. 2004) which has been further 
improved to Zt — 0.43 ± 0.07 (la) c.l. (Riess et al. 2007). Many authors have used these 
values of the transition redshift for imposing constraints on the cosmic parameters as an 
independent and trustworth discriminator for cosmology (Gardner 2005; Virey et al. 2005; 
Qiang & Zhang 2006; Gong & Wang 2006; Gong & Wang 2007, Xu et al. 2007; Ishida et al. 
2008). 

In this article, wc constrain the transition redshift by expanding the deceleration pa- 
rameter around its present value. A detailed analysis is carried out for 2 different parame- 
tcrizations: q = Qo + QiZ and q = qo + qiz{l + z)^^. Wc show that the kinematic analysis 
based on the SNLS (Astier et al. 2006), as well as the one recently compiled by Davis et 
al. 2007, yield a transition redshift Zt ~ 0.6 (lex) in better agreement with the fiat ACDM 
model {zt = {2VLf^/VLrnY/^ — 1 ~ 0.66) than the one provided by Riess et al. (2007) using the 
High-z Supernova Search (HZSNS) sample. 



To begin with, let us assume that the Universe is spatially fiat as motivated by inflation 
and the WMAP results (Spergel et al. 2007, Komatsu et al. 2008). Following Turner & 
Riess (2002) and Riess et al. (2004), the luminosity distance is defined by the following 
integral expression (with c—1) 



{z > 0.5). 



2. 



KINEMATIC APPROACH AND ANALYSIS 




(1) 
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Fig. 1. — Residual magnitude versus redshift is shown for SNe type la. The SNe of Riess et 
al. 2007 are represented by red circles, and, the remaining ones (Astier et al. 2006, SNLS) 
by open diamonds. Data and kinematic models of the expansion history are shown relative 
to an eternally coasting model, q{z) = 0. Models representing specific kinematic scenarios 
are illustrated. 

where the H{z) = a/a is the Hubble parameter, and, q{z), the deceleration parameter, is 
defined by 



Although generalizable for non-zero curvature, Eq. (1) is not a crude approximation 
as one may think at first sight. In the framework of a flat FRW type universe, it is an 
exact expression for the luminosity distance which depends on the expressions of the present 
Hubble constant, Hq, and the epoch-dependent deceleration parameter, q{z). The simplest 
way to work with the coupled definitions (1) and (2) as a kinematic model for the SN type 
la data is by adopting a parametric representations for q{z). As one may check, in the case 
of a linear two-parameter expansion for q{z) = qo + zqi (Riess et al. 2004), the integral (1) 
can be represented in terms of a special function as (see Appendix A) 




(2) 









(3) 
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where qo = q{z = 0) is the present value of the deceleration parameter, qi is the derivative in 
the redshift evaluated at z = 0, and 7 is the incomplete gamma function (see, Abramowitz & 
Stegun 1972) with the condition gi — go > must be satisfied (more details in the Appendix). 
By using the above expressions we may get information about qo, qi and, therefore, about the 
global behavior of q{z). Note also that a positive transition redshift, Zt, may be obtained only 
for positive signs of qi (the variation rate of qo) since qo is negative and the dynamic transition 
(from decelerating to accelerating) happens at q{zt) = 0, or equivalently, Zt = —qo/qi- 
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Fig. 2. — a) Contours in the qo — qi plane from 182 SNe of Riess et al. (2007). The two- 
parameter model of the expansion history is q{z) = qo + qiz, and the contours correspond 
to 68.3%, 95.4% and 99.7% c.l. as indicated. The best fit to the pair of free parameters is 
{loili) = (—0.7,1.62). b) Likelihood function for the past transition redshift. The best fit 
is., = 0.43^n^ (la) (2a). 



Another parametrization of considerable interest is q{z) = qo + qiz/{l + z) (Xu et al. 
2007). It has the advantage to be well behaved at high redshift while the linear approach 
diverges at the distant past. In this case we can write the integral (1) as (see Appendix A) 

Dl{z) = ^^^e'^^q;^'^^'^\^{q, + qo,q^) 

-i{qi + qo,qi/{i + z))], (4) 

where qo = q{z = 0) is the present value of the deceleration parameter, qi is the parameter 
yielding the total correction in the distant past {z ^ 0,q{z) = qo + qi)- Note also that a 
positive transition redshift, zt, may be obtained only for positive signs of gi > \qo\ since 
go is negative and the dynamic transition (from decelerating to accelerating) happens at 
q{zt) = 0, or equivalently, Zt = -qo/{qo + Qij- 
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The so-called Gold07 sample from the HZSNS team is a selection of 182 SNe la events 
distributed over the redshift interval 0.001 < z < 1.8, and constitutes the compilation of 
the best observations made so far by them which were completed by 16 events from the 
Supernova Cosmology Project observed with the Hubble Space Telescope (HST). 

The current data from SNLS collaboration correspond to the first year results of its 
planned five year survey. The total sample includes 71 high-a; SNe la in the redshift range 
0.2 < z < 1 plus 44 low-z SNe la. This data set is arguably (due to multi-band, rolling search 
technique and careful calibration) the best high-^ SNe la compilation to date, as indicated 
by the very tight scatter around the best fit in the Hubble diagram and a careful estimate 
of systematic uncertainties. Another important aspect to be emphasized on the SNLS data 
is that they seem to be in a better agreement with WMAP 3-years results (Spergel et al. 
2007) than the previous Gold04 sample observed by Riess et al. (2004). For a more detailed 
discussion see e.g., Jassal et al. 2006. 

For completeness, we also consider the Type la Supernovae compilation by Davis et al. 
2007. This extended sample is formed by 192 SNIa consisting of 45 SNe from a nearby SNe 
la subsample, plus 57 from SNLS and 60 intermediate redshift SNe from ESSENCE (Wood- 
Vasey et al. 2007), and, finally, 30 high redshift Gold07 SNe with internally consistent 
magnitude offsets. The supernovae span the redshift range 0.2 < 2; < 1.8 (Riess et al. 2007). 
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Fig. 3. — a) Contours in the qo — qi plane for 182 SNe type la data (a) Contours in the qo — qi 
plane for 182 SNe type la data {Gold07 sample) considering q{z) = qo + qiz/{l + z). Now, 
the best fit to the pair of free parameters is (go, qi) = (—0.84, 3). b) Likelihood function for 
the past transition redshift. The best fit is zt = 0.39lo:o6 (1^) laog (2a). 

In figure 1 we show the theoretical predictions of the kinematic approach to the residual 
Hubble diagram with respect to an eternally coasting Universe model {q{z) = 0). Such a 
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diagram is usually presented in terms of the bolometric distance modulus which is defined 

hy m — M = blogD^^z) + 25 (Peebles 1993). The different models arc characterized by the 
selected values of qo and gi, as depicted in the diagram. The SNe type la data shown in the 
panel comprise both the SNLS (Astier et al. 2006) and HZSNS (Ricss et al. 2007) samples 
as indicated. At this point, it is natural to ask about the likelihood contours in the plane 
(go,?i) and the probability of the transition redshift derived for each sample separately, as 
well as to the whole set of data. 

In order to obtain that, let us consider the maximum likelihood that can be determined 
from a statistics for a given set of parameters p 



where cTjuq. is the uncertainty in the individual distance moduli, is the dispersion in 
SNe redshifts due to pecuhar velocities (Riess et al. 2004, 2007), and the complete set of 
parameters is p = {Ho, qo, qi)- 

In what follows we investigate the bounds arising on the empirical q{z) parameters and 
the probability of the redshift transition for each SNe type la sample. By marginalizing the 
likelihood function over the nuisance parameter, Hq, the contours and the probabilities of 
the transition redshift for each sample are readily computed. 



In figure 2a we show the plane qo — qi for confidence levels of 68.3%, 95.4% and 99.7% 
for the Gold07 sample. Note that it strongly favors a Universe with recent acceleration 
{qo < 0) and previous deceleration {dq/dz > 0). With one free parameter the confidence 
region is —1.06 < qo < —0.42 and 0.6 < qi < 2.8 with (95.4%) confidence level. It should be 
remarked the presence of a forbidden region forming a trapezium. The horizontal line in the 
top is defined by gi = which leads to an infinite (positive or negative) transition redshift 
while the segment at 45° is the infinite future {zt = — 1). Actually, since z — Uo/a — 1, in 
the infinite future (a — > oo) one finds that 2; — > — 1. In addition, the values of Zt associated 
with the horizontal segment in the bottom are always smaller than -1 (therefore, after the 
infinite future!) since —1.5 < Zf < —1. Finally, one may conclude that the vertical segment 
is associated with Zt < —1.5, thereby demonstrating that the hachured trapezium is actually 
a physically forbidden region. 

In Figure 2b, one may see the probabihty of the associated transition redshift Zt defined 
by Q{zt) — 0. It has been derived by summing the probability density in the versus the 




(5) 



2.1. Gold Sample Analysis: q = qo + Qiz 
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dq/dz plane along lines of constant transition redshift, Zt — —qo/{dq/dz). The resulting 
analysis yields Zt = 0.431°;°^ (1 c" ) -ao9 (2*^ ) ^'^^ '^^'^ hcc parameter which is in reasonable 
agreement with the value Zt = 0.46 ±0.13 quoted by Riess et al. (2004). In our analysis, the 
asymmetry in the probability of Zt is produced by a partially parabolic curve obtained when 
is minimized. It is clear that the central value {zt = 0.43) does not agree with several 
dynamical fiat models in 68.3% and 94.5% confidence level. In particular, zt = 0.75 in the 
flat ACDM concordance model (fim = 0.27). 

2.2. Gold Sample: q — qo + qiz/1 + z 

In figure 3a we show the plane Qo — qi for confidence levels of 68.3%, 95.4% and 99.7% 
using the Gold07 sample for the above phcnomcnological law. Again, wc sec that a Universe 
with recent acceleration {qo < 0) and previous deceleration {dq/dz > 0) is strongly favored. 
With one free parameter the allowed region is —1.25 < qo < —0.15 and — 1.5<gi<4.9 with 
(95.4%) confidence level. Note also the presence of a physically forbidden region forming a 
trapezium (cf Fig. 2a). As one may check, the horizontal line at the bottom is defined by 
qi — which leads to a transition redshift in the infinite future {zt = —1) while the 45° 
segment also defines an extreme line of transition redshifts. Note also that all transition 
redshifts associated with the vertical segment of the triangle (on the axis defined by qi) take 
values on the future, more precisely, on the interval [-1,-1/2]. 

In Figure 3b we display the probability of the transition redshift zt which is defined by 
q{zt) = 0. The resulting analysis yields Zt = 0.39l|^:J^ (la) tofg i^^) for 

one free parameter. 

This central value agrees with the flat concordance model only at 95.4% {2a). Note that the 
straight lines denoting Zf = 0.3 (transition redshift for a fiat ACDM with ~ Q\ ~ 0.5) 
and Zt = 0.9 (for a ACDM with flm — 0.2 and fl\ ~ 0.8) are outside of the 2a regions in 
both parametrizations (see Figs. 2a and 3a). 

2.3. Analysis from SNLS: q — qo + q\Z 

Let us now consider the Astier et al. 2006 data sample to constrain the plane qo — qi 
for the above linear phenomenological description. 

In figure 4a we show the plane qo — qi with confidence levels of 68.3%, 95.4% and 99.7% 
to the SNLS data set. With one free parameter the confidence region is —1.17 < go < —0.16 
and —1.2 < qi < 3.4 with (95.4%). In the panel we sec that Zf = 0.3 is only marginally 
compatible at 3a while Zt = 2.0 is well inside with reasonable confidence {2a). Probably, 
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Fig. 4. — a) Constraints on the deceleration parameter from the SNLS sample with law 
q{z) = qo + qiz. Data favor the recent acceleration (go < 0) and past deceleration (gi > 0). 
Confidence intervals are showed in the diagram, b) Likelihood function for the past transition 
redshift. Due to the absence of data at high redshifts, we see that peak of the likelihood 
function for the SNLS sample is still pronounced but the upper limit of the transition redshift 
is not very well defined (confront with Fig. 2b). 



this apparent conflict comes from the fact that there are few data with high dispersions in 
z > 0.8 (see commentary by Astier et al. in their paper). Interestingly, the forbidden region 
now crosses the confidence level curves (even at Icr), however, the SNLS data (different from 
the GoldOT sample data) are compatible with the existence of a transition redshifts in the 
future, that is, beyond the forbidden region (2; > — 1). In particular, such a result reinforces 
the interest to examine higher order corrections (beyond the first order expansion) within 
the kinematic approach. 

In figure 4b we show the likelihood function for the transition redshift. Our analysis 
furnishes the best fit Zt = OMtlfi (la) for one free parameter. The likelihood is not well 
behaved at high redshifts because its upper limit is not known. Actually, although supplying 
an excellent adjustment when combined with the WMAP 3-years and cosmic concordance 
model, the superior limit of zt is not tightly constrained by the present SNLS sample. This 
poor inference should be expected from the very beginning in virtue of the absence of data 
for 2; > 1 in the SNLS sample. 
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Fig. 5. — a) Constraints on the deceleration parameter derived from the SNLS sample for 
q{z) = qo + Qiz/il + z). Again, the data favor the recent acceleration and past deceleration, 
b) Likelihood function for the past transition redshift. Note that the peak of the likelihood 
function is still pronounced but the upper limit of the transition redshift is not very well 
defined (confront with Figs. 3b and 4b). 

2.4. SNLS: q = Qo + Qiz/l + z 

In figure 5a we show the plane go — qi to the SNLS data set in the context of the above 
quoted parameterization. The allowed regions at 68%, 95% and 99% are indicated. With 
one free parameter the confidence region is —1.24 < Qo ^ —0.16 and —1.5 < qi < 4.9 with 
(95.4%). In the panel we see that Zt = 0.3 (fiat ACDM with Q^n - 0.5) is outside of the 
allowed {2a) region while Zt = 2.0 (fiat ACDM with fi^ — 0.07 and ~ 0.93) is well inside. 

It is worth noticing that the forbidden region is now represented by a right-angle triangle. 
Interestingly, it crosses the confidence level curves (even at la), however, the SNLS data 
(different from the Gold07 sample data) are compatible with the existence of a transition 
redshifts in the future, that is, for negative value of z. In particular, such a result reinforces 
the interest to examine higher order corrections (beyond the first order expansion) within 
the kinematic approach. In particular, this means that some theoretical problems associated 
with the existence of horizons in the future of models that accelerates forever must be solved 
in a natural way (see, for instance, Carvalho et al. (2006)). 

In figure 5b we show the likelihood function of the transition redshift. Our analysis 
furnishes the best fit Zt = 0.74 with inferior value 0.42 (lex) for one free parameter and the 
upper value to zt is too high. Actually, it is very high and completely incompatible with 
observations (see also Fig. 4b). Such a poor inference at high z from SNLS sample (Astier 



et al. 2006) for both parameterizations should be expected because the absence of data for 
z > 1. This is calhng our attention to the importance of observing high redshift supernovae 
within the SNLS collaboration. 




Fig. 6. — a) Constraints for the deceleration parameter are derived from the Davis et al. 
2007 analysis for q{z) = qo + qiz. 68.3%, 95.4% and 99.7% confidence intervals are showed in 
the graph, b) Likelihood function for the past transition redshift for a two-parameter model 
of the expansion history. Our analysis furnishes the best fit Zt = 0.60lo'^^ (Icr)- 



2.5. Davis 2007 Sample: q = + q^z 

Let us now discuss the constraints within the linear approach for the Davis et al. data 
set. In figure 6a we display the plane qo — qi for confidence levels of 68.3%, 95.4% and 99.7% 
as indicated. Note that such a sample also favors a Universe with recent acceleration {qo < 0) 
and previous deceleration {dq/dz > 0). With one free parameter the confidence region is 
— 1.0 < qo ^ —0.36 and 0.11 < qi < 2.2 with (95.4%). In the same panel we also see that 
Zt = 0.3 is completely outside even of the significant 3a region while Zt = 0.9 is well inside 
with high confidence. 

In figure 6b we show the likelihood function for the transition redshift. We can see 
that Zt = 0.60^011 (la) for one free parameter. These results are in better accord with flat 
LCDM models. Such a fact may be somewhat related to the higher degree of statistical 
completeness obtained when the SNe type la samples are added. 
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Fig. 7. — a) Constraints on the deceleration parameter are derived from the Davis et al. 2007 
analysis. The data favor the recent acceleration (go < 0) and past deceleration (gi > 0) with 
high confidence for q{z) = Qo + qiz/i). + z). b) Likelihood function for the past transition 
redshift for a two-parameter model of the expansion history, q{z) = qo + qiz/{l + z), from 
SNe la. Our analysis furnishes an accurate best fit Zt = 0.38^o;Qg (Icr). 



2.6. Davis 2007 Sample: q = qo + qiz/1 + z 

In Figure 7a we display the plane qo — qi when Davis sample is considered for q{z) = 
qo + qiz/{l + z) parameterization. With one free parameter the confidence region is —1.17 < 
qo < —0.36 and 0.39 < gi < 3.99 with (95.4%). In the same panel we also see that zt = 0.3 
is outside while zt = 0.9 is well inside with high confidence. 

In figure 7b we show the likelihood function of the transition redshift. We can see that 
Zt = 0.38lo Qg (l(j) for one free parameter. This results is comparable with the constraints 
from the GoldOl sample for the same parameterization (see subsection 2.2). It is also inter- 
esting that the Davis et al. sample also suggests (for both parameterizations) the existence 
of a transition redshift in the future {z < 0). Such a possibility it will be discussed with 
more detail in a forthcoming communication (Lima and Cunha 2008). 



3. RECONSTRUCTING THE DECELERATING PARAMETER 

At this point it is interesting to investigate how a joint analysis constrains the redshift 
evolution of the decelerating parameter. However, since the time evolutions are quite similar 
for the different parameterizations, in what follows we focus our attention on the Gold07 
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Fig. 8. — a) Evolution of the decelerating parameter as a function of the redshift based on 
the expansion q{z) = qo + Qiz. b) The same of Fig. (a) for a two-parameter expansion 
history driven by q{z) = Qo + Qiz/^l + z). In the two panels the shadowed regions means 
2a level for the Gold07 sample. The dotted horizontal lines represent the coasting model 
{q{z) =0). In both parameterizations the data favor a low redshift accelerating stage of the 
Universe (go < 0) and a past decelerating phase (gi > 0) with high confidence level. 

sample. The basic results are displayed in Figures 8a and 8b. 

In Figs. 8a we see the evolution of the decelerating parameter as a function of the 
redshift for the parameterization q{z) = qo + qiz. The shadow denotes the 2a region. 

In Fig 8b the evolution is shown for a two-parameter expansion history driven by q{z) = 
qo + qiz/{l + z). For both cases, we also show the lower and upper values {2a) of the 
transition redshifts which is defined by the condition {q{z) = 0). Note also that for both 
parameterizations the data favor the recent acceleration (go < 0) and past deceleration 
(gi > 0) with high confidence level. 

In Table 1 we summarize the basic constraints on the transition redshift {zt) established 
in the present work. As shown there, for a given phenomenological law the limits were 
derived separately for each sample of SNe type la data. 

4. CONCLUSION 

The most remarkable dynamical feature distinguishing dark energy and matter domi- 
nated models is the present accelerating stage. However, the observed structures must be 
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formed in a past decelerated phase dominated by matter. This leads naturally to a dynamic 
phase transition from decelerating to accelerating in the course of the Universe's evolution. 
This is the most natural scenario within the framework of the general relativity. It has been 
confirmed by several disparate observations ranging from the anisotropy of the CMB, SNe 
la, galaxy clusters, age of the universe, however, even considering the high-quality of the 
present data the whole set is still compatible with many different dark energy models. 

The basic information about the deceleration parameter has been obtained here through 
a model independent analysis using two different phenomenological laws, namely: q{z) = 
qo + qiz, as proposed by Riess et al. (2004) and q{z) = qo + qiz/{l + z) (Xu et al. 2007). 
The contours to the plane (go, qi) and the likelihood function for transition redshift were 
discussed in detail with basis on the data sample by Riess et al. (2007), Astier et al. (2006) 
and Davis et al. (2007). In order to compare the predictions, we have presented the analysis 
first for each sample separately. As we have seen, although considering the presence of no 
transition regions (see Figs. 2a, 3a, 4a, 5a, 6a and 7a), both sample strongly favor a Universe 
with recent acceleration {qo < 0) and previous deceleration {dq/dz > 0). 

The likelihood function for the transition redshift is more well defined for the Riess et 
al. (2007) and Davis et al. (2007) data. In table 1 we summarize the main results related to 
the transition redshift Zt in the present kinematic approach. For both samples, our analysis 
provides a model independent evidence that a dynamic phase transition from decelerating 
to accelerating happened at redshifts smaller than unity. Hopefully, the constraints on Zt 
will be considerably improved in the near future with the increasing of supernovae data at 
intermediate and high redshifts. 
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Appendix 

A. Distance Luminosity in the Kinematic Description 

In this Appendix, we shall outline the method used to obtain the theoretical expression 
for Dl{z) when the kinematic approach is based on the expansion: q{z) = Qo + Qiz. The 
result is valid for a geometrically flat model. By definition: 



Ddz) = {l+z) f ^ (Al) 
Jo H{u) 

where we have set c = 1 and H is the Hubble parameter which in the kinematic approach 
must be obtained in terms of the decelerating parameter: 



ad dH ^iz) 



a dH 
H da 



By integrating the above equation and inserting the result into flAll) it follows that 

(1 + ^) 



Dr Z) 



Jo 



exp 



[1 + q{u)]d\n{l + u) 



du 



which is just ([T]). Replacing the expansion of q{z), we obtain 

Dl{z) = (1 + z)H-^ f (1+ uy-^ exp [-H du 
Jo 

where fJ^ = qi and u = qi — qo- Now, a simple variable change results 



Dr(z) 



e''qr"''[l{qi-qo,{l + z)q,: 



Ho 

-liQi - (lo,qi)] 



(A2) 
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where 7 is the incomplete gamma function (Abramowitz & Stegun 1972). 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



Interestingly, in the linear approximation the luminosity distance can also be directly 
written in terms of the transition redshift: 

-"0 V Zt 



7 



\qo\{l + l/z^)}-^{l + z) 



7 



(A6) 



For the second parametrization q{z) — qo -\- qiz / {1 -\- z) , we obtain 

Dl{z) = (1 + z)H-^ / (1 + uf-^ exp [-u^i] du 
Jo 



(A7) 



where ^ = qi and i' = qi + qo- As well as in the last integral, a simple variable change results 



Dl{z) = ^^e«^gr^'''^'^^[7(?i + ?o,?i) 
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-7(^1 + ?o,gi/(l + ^))]- 



(A8) 



Or yet, in terms of the transition redshift: 

(l + z) W2i(l±Hlf\qo\{l + Zt)y°^'' 
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Table 1. Limits on the transition redshift Zt 



Sample Best fit Confidence region {la) 



g = 9o + qiz 

RiessOT 0.43 0.38 < zt < 0.52 

AstierOe 0.61 0.40 < zt < 4.29 

DavisO? 0.60 0.49 < zj < 0.88 

<-! = <io + '71-/(1 + z) 

RicssO? 0.39 0.33 < zt < 0.49 

AstierOe 0.74 0.42 < Zt high 

Davis07 0.38 0.32 < Zt < 0.55 



